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Abstract 

The 1-parameter family of transformations identified by Barbero and Immirzi plays a significant role in 
non-perturbative approaches to quantum gravity, among them Loop Quantum Gravity and Spin Foams. It 
facilitates the loop quantization programme and subsequently the Barbero-Immirzi parameter (7) arises in 
both the spectra of geometrical operators and in the dynamics provided by Spin Foams. However, the debate 
continues as to whether quantum physics should be Barbero-Immirzi parameter dependent. Starting from a 
5^ , discrete SO(A)-BF theory phase space, we find two possible reductions with respect to a discrete form of the 

simplicity constraints. The first reduces to a phase space with 7-dependent symplectic structure and more 
generally in agreement with the phase space underlying Loop Quantum Gravity restricted to a single graph 
- a.k.a. Twisted Geometries. The second, fuller reduction leads to a 7-independent symplectic structure 
■ on the phase space of piecewise-flat-linear geometries - a.k.a Regge geometries. Thus, the 7-dependence of 

' physical predictions is related to the choice of phase space underlying the quantization. 
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; 1 Introduction 

OV 

Since its inception, the study of the gravitational force has proved to be a fruitful yet frequently perplexing 
CSJ ■ endeavour. With the advent of quantum theory and its manifold applications to the description of matter, the 

quantization of gravity seems a natural step towards a more cohesive picture of the world around us. Holding 
true to form, however, the twin concepts of dynamical geometry and diffcomorphism invariancc appear to lie 
at the root of gravity's stubborn resistance to a fundamental quantum mechanical description. 

For quite some time now, there have existed a number of approaches to this challenge that, in one manner 
or another, involve replacing the continuum manifold by a simplicial one. We shall mention just two of these: 
quantum Regge calculus [1, 2, 3] and the more modern Spin Foam models [4, 5]. The former is a Lagrangian 
path integral approach where all the paths have support on a fixed simplicial manifold and are weighted by 
the exponential of the Regge action - a discretization of the Einstein-Hilbert action. The latter are state sum 
models that can be interpreted in various ways (a property that makes them all the more appealing) and using 
different kind of variables, for instance non-commutative first order variables [6, 7], twisted geometry variables 
[8] or just holonomies [9, 10]. Spin foam models can also be viewed as a Lagrangian path integral approach 
[11, 12, 13], where once again all the paths have support on a fixed simplicial (or even cellular) manifold and are 
weighted by the exponential of a discretized form of the Hoist-modified Plcbanski action [14, 15]. The presence 
of the Hoist term not only entails that modern spin foam models have a dependence on the Barbero-Immirzi 
parameter (7) [16, 17, 18]. We shall present more details on this a bit later, but for the moment, let us note that 
quantum Regge calculus has no 7-dependence. More importantly, the ultimate significance of this parameter 
remains unclear. 
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Our aim in this paper is to ameliorate this state of affairs by clarifying the role of the Barbcro-Immirzi 
parameter in the classical discrete theory. 

But first, let us begin in the continuum, where a somewhat similar situation exists. Tersely speaking, Loop 
Quantum Gravity (LQG) is a traditional canonical quantization procedure applied to gravity that has provided 
a rigorous definition of the kinematical Hilbert space, along with an anomaly-free action of the Hamiltonian 
constraint operator [19, 20, 21]. These powerful successes stem from the ability to parameterize the phase space 
of continuum spatial geometries in terms of an SU(2) connection and a triad field [22], the same variables that 
parameterize the phase space of SU(2)-i3F theory [23]. Of course, these variable are still subject to further 
constraints: the rotational Gauss, 3-diffeomorphism and Hamiltonian constraints. The gauge-theoretic profile of 
this phase space is not a priori apparent, rather it is uncovered via a canonical transformation [16]. Remarkably, 
there is a 1-parameter family of such transformations, labelled by the Barbero-Immirzi parameter. 

While the classical theory is, at the end of the day, insensitive to these 7-dependent transformations, the 
Barbcro-Immirzi parameter subtly finds its way into the quantum theory [24] . The point is that transformations 
with different values of 7, define different pairs of conjugate (kinematical) observables. These are subsequently 
quantized, that is, they are represented as operators on a Hilbert space. The quantization procedures involves 
a choice of holonomics as basic variables, in this case, exponentiated (7-dcpcndent) connections. This leads to 
the Barbcro-Immirzi parameter manifesting itself in the Poisson brackets between basic classical variables and 
therefore in the commutators involving basic quantum operators. 1 

This quantization procedure even leads to the appearance of 7 in the Poisson and commutator brackets 
between spatial geometric quantities, the 3d-dihedral angles for example, despite 7 having been introduced only 
in the connection (which encodes the extrinsic curvature). In particular the spectra of geometrical operators, 
notable the area and volume operators [28, 29, 30], turn out to depend on 7 and thus show that quantum 
theories for different values of 7 arc unitarily incquivalcnt. 

This is not the end of the story in the continuum, however. One may view the path to this loop gravity phase 
space within the context of the Hoist- modified Einstein-Cartan formulation [14]. Therein one begins with an 
auxiliary phase space parameterized by bi- vectors and a Lorentz connection. These variables are identical those 
of a Lorentzian BF theory and are subject to the famous simplicity constraints (both primary and secondary) 
and Lorentz Gauss constraint, in addition to the 3-diffeomorphism and Hamiltonian constraints. 

Importantly, the primary and secondary simplicity constraints are 2nd class and therefore in any quantization 
procedure, they should be dealt with at the classical level. There are two completely equivalent methods to 
reduce with respect to these constraints: i) explicit solution or ii) Dirac bracket construction. 

To mention only the pertinent details, explicit solution yields a spatial Lorentz connection, whose construc- 
tion involves 7, and its conjugate momentum. Moreover, a subsequent gauge-fixing of the boost part of the 
Gauss constraint lands us back in the phase space underlying the LQG approach [31, 32]. Thus, one expects that 
the Barbcro-Immirzi parameter still plays a fundamental role after quantization. On the other hand, by con- 
structing a Dirac bracket, one can parameterize the reduced phase space using yet another Lorentz connection 
(and its conjugate momentum). Remarkably, the Barbero-Immirzi parameter has disappeared, i.e. it does not 
show up in the Poisson brackets between this connection and conjugated variables nor in the other constraints 
[33, 34, 35]. Therefore, any subsequent quantization would inherit this property. Given the presence of the 
Lorentz connection, it was hoped that one could perform an analogous loop quantization [36]. Unfortunately, 
this connection is non-commutative (w.r.t. the Dirac bracket) and as a result, only a formal quantization has 
been achieved. Despite that, it has spurred on much debate as to whether quantum gravitational physics is 
7-dependent [37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47]. 

In particular, let us apply the Jacobi identity, which holds for both Poisson and commutator brackets, to {e 79 , {pi,J>2}}; q,P are 
analogues of the connection and triad, respectively. One notices that a non-vanishing bracket between momentum variables implies 
{pi,P2} ~ 7, assuming that {q,Pi} = 1. Indeed, the necessity of non-commuting fluxes in Loop Quantum Gravity [25, 26, 27] 
follows from a similar kind of argument. 
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As mentioned earlier, the apparent dilemma with respect to 7-depcndcnce has a counterpart in the discrete 
setting. We shall reveal it along the following lines. 

Despite the progress in the LQG programme, the analysis of the Hamiltonian constraint operator in the full 
kincmatical Hilbcrt space has remained for the most part at a rather abstract level; one can formally define 
a projector onto the physical state space. As a result, attention has concentrated on providing a definition of 
the projector onto physical states for spin networks with support on a fixed graph. In essence, this corresponds 
to a truncation of the theory, with some very interesting repercussions. One considers a spin-network as the 
dual graph to a cellular decomposition of the spatial 3-manifold. Furthermore, there is a classical phase space 
underlying that part of the Hilbert space associated to each spin network graph. As one might imagine from our 
earlier discussion, it is essentially the phase space of discrete S\J(2)-BF theory. However, it may be recast in a 
more geometrical fashion known as Twisted geometries [48] and is a discretization of the phase space underlying 
LQG. Thus, as we shall detail later on, the Poisson brackets between basic variables, and the constraints 
have explicit Barbero-Immirzi parameter dependence, leading again to the appearance of 7 in the spectrum 
of geometric operators. Consequently, the standard spin foam quantization procedure results in 7-depcndent 
amplitudes weighting the paths. 

However, there is another method to quantize gravity in the discrete setting known as quantum Reggc 
calculus, which is 7- independent. Recently, a phase space has been prescribed for the classical theory [49, 50]. 
This phase space describes Regge geometries geometries or more precisely, piccewisc-lincar-flat geometries. 

It should be noted that Twisted geometries are not identical to Regge geometries. Regge geometries are 
a more traditional interpretation of discrete geometry in that one has an unambiguous assignment of edges 
lengths and deficit angles to the simplicial 3-manifold. This is not case for Twisted geometries; while one can 
construct edge lengths and deficit angles, they are dependent on a choice of tetrahedron and edge respectively, 
whose data are used to construct these geometrical quantities. 

As a consequence, an important statement can be made at this point. In the continuum, the variables of 
SV(2)-BF theory capture spatial geometries, as well as the conjugate extrinsic curvatures. In the discrete, they 
capture a generalized class of spatial discrete geometries (and conjugate variables). The discrete SU(2)-BF 
phase space, as well as the continuum one, lead to 7-dcpendcnt Poisson brackets between basic variables and 
consequently to 7-dependent spectra of geometric quantum obscrvablcs. This opens the question to whether 
a further reduction of the discrete phase space to a phase space describing Regge geometries is possible and if 
so, whether such a description would involve 7 in any way. If the latter were not the case, then 7 could not be 
involved in the spectra of any geometrical observables. 

We shall see in this paper how the different discrete geometrical phase spaces for are related to each other. 
More precisely, this paper summarizes relevant findings of two previous papers [51, 52] on the various forms of 
discrete simplicity constraints and gives the analysis of the various Dirac bracket structures, so as to infer the 
status of 7 in the different phase spaces. In Section 2, we give a very brief description of the continuum Hoist 
phase space. Thereafter, in Section 3, we discrctize the spatial manifold along with the dynamical system. In 
this canonical setting, this means that we provide a discrete counterpart of the phase space variables, along 
with an appropriate symplectic structure subject to a discrete version of the constraints. As we shall deal 
with Ricmannian gravity, we shall have at the outset an SO(4)-£?F theory phase space subject to a discrete 
form of the Gauss, simplicity and diffeomorphism constraints. Our analysis, in Section 4, then concentrates on 
completing a reduction with respect to the Gauss and simplicity constraints. 

The upshots of this analysis are the following: 

- If we reduce by a subset of the discrete simplicity constraints, we arrive at the phase space of Twisted 
geometries, in its 7-dcpcndcnt description. 

- If we reduce by all the discrete simplicity constraints, we arrive at the phase space of Regge geometries, 
with 7-independent Poisson brackets between basic variables. 
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We can draw a number of conclusions: 

- In effect, the above results highlight a very common property of the discretization mechanism: at least at 
the outset, one has a certain amount of freedom when it comes to discretizing a dynamical system. 

Here, one such choice is to hrst reduce the Hoist phase space in the continuum down to the loop gravity 
phase space and then discretize. In this fashion, one is naturally led to the phase space of Twisted 
geometries. 

Another choice is to first discretize the larger Hoist phase space, based on 50(4) holonomies. In this case, 
one has some freedom in how one discretizes the simplicity constraints. Of course at this kinematical 
stage, such choices are based on prejudice. If one is motivated to maintain strict discrete geometricity, 
one chooses a larger set of discrete constraints and one arrives at the phase space of Regge geometries. 
However, if one is motivated to get to twisted geometries, then one chooses a smaller set of the discrete 
constraints and one finds 7 dependent Poisson brackets between basic variables. Thus the appearance 
of 7 can be seen due to an enlargement of the discrete configurations, on which a dynamics is defined. 
Ultimately, the question will be whether one can find a consistent dynamics based on these configurations, 
leading to the correct classical limit, in which 7 should not have any physical significance. 

- The respective appearance and disappearance of 7 in the Twisted and Regge phase spaces and their cor- 
responding quantum theories might not be too surprising. The constraints which would reduce the phase 
space, equivalent to Twisted geometries, to the Regge phase space are second class. Following standard 
Dirac procedure, one should solve these constraints classically and perform quantization afterwards. 2 

- Finally, our analysis highlights a basic difference between the 7-dependence argument in the continuum 
and this one in the discrete. In the continuum, the 7-dependent description of the reduced phase space 
is classically equivalent the 7-independent one. They are two parameterizations of the same phase space, 
related by a canonical transformation, albeit a highly non-trivial one. In the discrete, the 7-indcpcndent 
Regge geometry phase space is truly a subspace of the 7-dependent Twisted geometry one. Furthermore, 
all the 4d-dihedral angles commute with respect to the Dirac bracket, while their continuum counterpart 
in the 7-independent description, the Lorentz connection, is non-commutative. 



2 Plebanski's action and phase space 

The Plebanski action for gravity on a smooth manifold Ai takes the form: 



S 



Pleb,M 



[X, to, . 



d A x 



M 



M VP < ii) Xxil ,F[w} vp > +i ^ < X x ^X vp y 



(1) 



We have introduced co-ordinate charts on space-time parameterized by x^ = (t,x l ), where fi G {0, 1, 2, 3} and 
i € {1, 2, 3}. The dynamical fields are an so(4)-valued bi- vector field X and an so(4)-connection to. < , > and 
-< , >- arc the two bilinear forms on so(4). 3 Moreover, the action contains a map : so(4) — > so(4) = S + 



2 Yet another interesting alternative is [53, 54, 55], which shares with this approach an SO(4) phase space. But instead of 
implementing primary and secondary simplicity constraints, only the primary are implemented. For a consistent dynamics one 
needs to adjust the Hamiltonian constraints. For a discussion of the various ways and difficulties to implement primary simplicity 
constraints (which turn out to be non-commuting in the quantum case despite being first class classically) see [56] . 

3 The curvature of this connection is F[w] pv A = di fl w 1/ -i A +C A bc w [^ F 
Meanwhile, the bilinear forms arc defined: 



< X ,Y >= X A S AB Y 1 



and 



c where C ABC are the structure constants of so (4). 

(2) 



< X, Y y = X A e AB Y B 
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where 7 e 1 is the Barbcro-Immirzi parameter. Thus, the first term in the action describes a topological theory 
known as BF theory (with Hoist modification [14]). The Lagrange multiplier (f> x vP is responsible for imposing 
the Lagrangian simplicity constraints. 

As we outlined in the introduction, to perform a phase space path integral quantization of the above system, 
one needs more information about the constraint structure of the above theory. In fact, a detailed analysis of 
the canonical structure has been performed [57, 58], including an explicit symplectic reduction. One finds that 
at the outset, the constraint structure is rather intricate, but after a reduction with respect to a subset of the 
constraints, one arrives at a more familiar phase space Vy, (where £ is the characteristic leaf of the space-time 
foliation: M = lx S). Its parameterization and symplectic structure are: 

(V(if), ™W B (y) = ^ 5bc h) X jk c {y)) and { Wl A (x), ™TP B (y)} = V 5 A B 6 {3) (x,y) , (3) 

Thus, at this stage of the proceedings, we have the phase space of SO(4)-£?F theory. However, the phase space 
of gravity is a hypcrsurface within this phase space. In other words, these phase space parameters are still 
subject to the following constraint system: 



Gauss: 


Q A 




Primary simplicity: 




= ^n^iFV, 


Secondary simplicity: 




= e " m < 4 < h mp WlF, A.IP') > , 


Vector: 


K 


= e 0ijk W A S AB F[w] k B , 


Scalar: 


Ao 


= hijlPAC^ F[wYb ■ 



We have split the constraints here. The vector and scalar constraints are related to the imposition of 4- 
dimensional diffcomorphism symmetry, while the Gauss, primary and secondary simplicity constraints reduce to 
a 3-dimcnsional geometrical phase space. We shall impose the Gauss and simplicity constraints in the discrete 
setting, while leaving the dynamical constraints for future consideration. (In fact, in previous work [51], we 
have examined the imposition of a more tractable "flat dynamics" constraint.) 

Remark: For completeness, we note here that this phase space is essentially the same as that of the Holst- 
modihed Einstein-Cartan action [58]. Moreover, one arrives at the phase space of LQG in the following way: 
one gauge-fixes the boost part of the above Lorentz Gauss constraint, which allows one to explicitly solve all 
the simplicity constraints in terms of an su(2)-connection A and a triad field E. The parameters and symplectic 
structure on this phase space are: 

(A 4 %t) ee T*[E{x)] - jKi a (x), Eh(y) = TPabdri) and {^"(5), E\{y)} = ySJ S\ 6^(x, y) . (5) 

where T denotes the spatial Levi-Civita connection constructed from the triad and K encodes the extrinsic 
curvature tensor. As one can see, this phase space is essentially identical to that of SU(2)-BF theory, although 
it is still subject to the vector and scalar constraints. Moreover, it is the discretization of this phase space that 
leads to twisted geometries. 

4 hij is the spatial metric defined by: hh lJ = < IP, IP >, h = dct(hij) where h tJ is the inverse of hij. 
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3 Discrete Plebanski phase space: Va 3 



As detailed comprehensively in [51, 52], the discretization procedure S — >• A3 entails the replacement of Vy, 
with Va 3 along with their corresponding constraint sets. On the simplicial complex, tetrahedra arc indexed by 
{i}. Therefore, triangles (edges) are indexed with respect to the couple (a triple) of tetrahedra to which they 
belong. We illustrate these choices in Figure 1. 




Figure 1: Indexing the various sub-simplices of A3. 



The variables parameterizing Va 3 arc a discrete connection and a discrete bi- vector field, denoted by (M^-, X^), 
respectively. We illustrate this in Figure 2. 





Figure 2: The discrete bi- vector and discrete connection. 



The symplectic structure on V a is: 



{Mij± AD , X l]± B } 
{X ij ± A ,X ij± B } 



{M ij± AB ,M ij± UD } 



7± 1 



C AB C M ij± CD 



1 



7± 1 




r<AB y C 

k c Xij± 



(6) 



where we have separated the degrees of freedom into the self-dual and anti-self-dual sectors using the appropriate 
projectors on so (4). More details can be found in Appendix A. 

At this moment, we have an important choice to make concerning the discrete constraint. To phrase our 
dilemma differently, what properties should we like the discrete reduced phase space to possess. Naturally, 
we would like to arrive ultimately at a discrete geometrical phase space and we take the view that the most 
appropriate choice is the set of piecewise linear flat geometries - known also as Rcggc geometries in the literature. 
We shall take this as our guiding principle when choosing a constraint set. 

With this in mind, there are several constraints to implement on this system if one is to remove non- 
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geometrical configurations. 



They take the form: 



Matching: 

Closure: 
Diagonal simplicity 
Cross simplicity: 
Edge simplicity: 



1 M lj ± A 
2 M lJ± AB 



Xji± + Mji± b Xij± 
(M ij± M Jl± ) AB - 6 AB 



Aij + 
<t>w+ - 

9ij;kl+ 



A 



ij- 
kjk- 

Oij-ki- 



kill the doubling of variables 
associated to each triangle, 



(7) 

The simplicity constraints are written in terms of gauge invariant quantities, detailed in appendix B, which 
acquire a geometrical meaning on the hypersurface specified by the above simplicity constraints. Aij denotes 
the area of the triangle {ij}, cos tfiij^ corresponds to the 3d-dihedral angle between the faces {ij} and {jk}, 
while cos8ij;ki is the 4d-dihcdral angle between the tetrahedra {i} and {j}. We provide an illustration of these 
quantities for geometrical configurations in Figure 3. An important point for the following is that there are three 
edge simplicity constraints and, a priori, three dihedral angles per triangle (and chiral sector). This abundance 
of different ways to define dihedral angles is indicated in the index structure Oij-u , where ij specifies the triangle 
in question and kl the particular edge which is used for the definition of the dihedral angle. Only on the final 
reduced (Regge) phase space these different dihedral angles per triangle are forced to agree. 







{i :l k} = {jt.l} 



Figure 3: The gauge invariant variables. 



Remark: We proved in [51] that this constraint set is sufficient. Namely, upon imposing the simplicity 
constraints, one arrives at the phase space, in which one can consistently reconstruct edge vectors. Furthermore, 
these edge vectors are properly parallel-transported by the discrete connection. Thus, in the gauge-invariant 
subspace, one can reconstruct a unique set of edge-lengths and their conjugate deficit angles; the reduced phase 
space is the phase space of Regge geometries. 

Remark: In [52], we proved that this constraint set, in gauge invariant form (also known as quadratic form of 
the simplicity constraints), is equivalent to a gauge variant or linear form of the simplicity constraints. In time 
gauge, this linear form of the simplicity constraints resembles quite obviously a discretization of the primary 
and secondary simplicity constraints of the continuum, equivalent to the reality conditions for the complex 
(self-dual) Ashtekar variables, i.e. E + = E_ and A + + A- = T. Here, E denotes the triad variables, and A the 
connection variables in the different duality sectors. T is the spatial Levi-Civita connection. 
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3.1 Discrete constraint analysis 

One must be careful when referring to the ^-angles as 4d-dihedral angles. A 4d-dihedral angle, 6ij is the angle 
between the normals to the tetrahedra {i} and {j}. Thus, it is labelled by the shared triangle {ij}. Off the 
constraint hypersurface, however, there are three different ways to define this angle (in terms of the discrete 
bi- vectors and connection). Hence, we use the nomenclature 9ij-ki which denotes the dependence on the edge 
{ijk} = {jil}- It is only on the geometrical phase space that these three definitions coincide. 

Pursuant to this remark, we shall pay some special attention to the simplicity constraints. It turns out that 
they have a rather subtle structure. At the outset, they possess a powerful yet simple appearance imposing: 

"LEFT-HANDED GEOMETRY = RIGHT-HANDED GEOMETRY" 

In fact, they also implicitly ensure that conditions are satisfied among quantities within one and the same 
chiral sector. To begin, one can show that there are tidy relationships among the 0-anglcs associated to a given 
triangle {ij}: 

@ij;k'l'± — 6ij;M± = &ij;kk'± — a ji-l'l± ■ (8) 

The quantities O^-m and 9ij-,k'i' ar e two definitions of the 4d-dihcdral angle associated to the triangle {ij}, 
corresponding to two different edges. In addition, oiij-kk' an d otji,w are the two definitions of the same 2d- 
dihcdral angle in the triangle {ij}, as computed within the two different tetrahedra sharing this triangle. 
Upon imposing the simplicity constraints, one finds that: 

4>ijk+ — 4>ijk- 

9ij;kl+ = —Oij;kl- } — ^7— > %;fe'i'+ = %:fe/+ and aij :k k>+ = j (9) 

cq. (8) 

&ij;k'l'+ = —8ij;k't>- J 

which means that on the constraint hypersurface, the 0-angles are really 4d-dihedral angles. They only depend 
on the triangle {ij}. Furthermore, equation (9) shows that if it should suit our purposes, we may replace some 
of the edge simplicity constraints with constraints among the a-angles in one chiral sector. They are known as 

the GLUING CONSTRAINTS [59, 51]. 

4 Imposing the constraints 

To begin, let us calculate the dimension of the initial phase space. For a generic triangulation A, the dimension 
of Va is 24Nf, where N t , Nf and N e are the number of tetrahedra, triangles and edges, respectively. 

Our next task is to find an irreducible set of constraints. Given that the reduced phase space is parameterized 
by Rcgge geometries, we know that this irreducible set must number: 24JV/ — 2N e . 
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4fc OF CONSTRAINTS 


PHASE SPACE (DIMENSION) AFTER REDUCTION 


Stage 1: 


1 M ij±1 2 Mjj± 


127V/ 


l2N f 


Stage 2: 




6A^ 4 (= 3N f ) 

v \ J / 


6Nf Gauge-invariant phase space 




Dji 
L j 


N f 




Stage 3: 




Nf 


3A/ Twisted geometries 




Eik.Jl 


N f 






Fijfc 


2N t (= Nf) 




Stage 4: 


A 


N f ~N e 


2N e Regge geometries 




T 


N f -N e 





(10) 



As one can see, we have partitioned the constraints into four separate subsets: STAGE 1 involves the matching 
conditions 1 Mij± and 2 Mij±; STAGE 2 involves the closure constraint Gi±; STAGE 3 involves the diagonal, cross 
and a subset of the edge simplicity constraints; while STAGE 4 involves the rest of the edge simplicity constraints, 
albeit in an alternative form. We shall detail their exact form later in this section. 

Remark: As one might imagine, we re-express the final edge simplicity constraints in another form in order 
to facilitate their imposition. Note that the gauge-invariant phase space possesses a basis with 6iV/ elements 
split half-and-half between the chiral sectors. The irreducible set of simplicity constraints contains 6Nf — 2N e 
elements. Thus, while we can impose 3Nf of these in their "left-handed = right-handed" form in STAGE 3, 
there is a clear obstruction to the imposition of the final 37V/ — 2N e in this manner. For this reason, we rewrite 
them as constraints on a single chiral sector, that is, as gluing constraints. 



4.1 Stage 1: Matching constraints 

The constraints are 2nd class and couple the {ij} variables with {ji} variables. As a result, we may deal with 
them in two equivalent ways: either we solve the constraints explicitly within the Poisson bracket, or we perform 
a symplectic reduction to obtain a Dirac bracket. We choose the former here as the constraints are easily solved: 



{M ij± AD ,X ij± B } 
{Xij± A , Xij± B } 



7 



7± 1 



7± 1 



C ABC M ij± CD 
r'ABC v , c 



{M ij± DA ,X ji± B } 



7 



{M ij± AB ,M ij± CD } = 0, 



7± 1 

{X%j± ,Xji± B } = 0, 
{M l0± AB ,M ]%± CD } = 0. 



C ABC M 13± DC , 



The number of matching constraints is 12 Nf so that the reduced phase space has dimension 12 Nf. 



(11) 
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4.2 Stage 2: Closure constraint 



The closure constraints Gy. are 1st class and generate SO (4) gauge transformations. We can take these con- 
straints into account by restricting to the constraint hypersurface and parameterizing the gauge orbits using a 
basis of gauge-invariant variables. The number of constraints is 6N t = 3Nf but since they are 1st class, they 
remove twice as many phase space dimensions. Thus, the resulting gauge-invariant phase space has dimension 
6Nf. To parameterize it, we shall use the following gauge- invariant parameters: 





# IN A3 


# USED 


Areas: 


Aij± 


2N f 


2N f 


3D-DIHEDRAL ANGLES: 


<f>ijk± 


12N t 


4N t = 2N f 


Averaged 4d-dihedral angles: 


6ij± 


2N f 


2N f 



(12) 



In each chiral sector, we use all the areas, two (non-opposite) 3d-dihcdral angles from each tetrahedron, and 
one averaged 4d-dihedral angle per face. It is a matter of taste that prompts us to use the average of the three 
4d-dihcdral angles attached to the triangles. We would rather not single out one of the three as being special. 
The rest of the gauge-invariant quantities can be rewritten in terms of these chosen parameters, so we have a 
basis (see Appendix B.2 for details). 

To keep matters clear, let us present just a schematic of the symplectic structure on the gauge-invariant 
phase space: 



{A±,A ± } 



{A±,0±} = ±2 



7± 1 



{A±,<f>±} 



W±,0±} = ±2 



7 /(i)({0±l) 
7± 1 A± 



{o±,e±} 



± 2- 



= ±2 



7 sma± 
-y±l A± 

7 /(2)({0±» 

- 7±1 A± 



where we present a much more comprehensive list in Appendix B. 



4.3 Stage 3: "left-handed = right-handed" simplicity constraints 

The simplicity constraints Dy, Cijk and Eik,ji are 2nd class. As we detailed earlier, we split them into two 
subsets. We impose the first subset at this stage. They enforce that the left-handed basis elements match up 
with their right-handed counterparts in the following fashion: 







# IN A3 


# USED 


Diagonal simplicity: 


Dij -- 


- Aij+ Aij- 


Nf 


N f 


Cross simplicity: 


Gijk = 


- 4>ijk+ — 4>ijk- 


6N t 


2N t = N f 


Partial edge simplicity: 




= e ij+ + 9ij- 


3N f 


N f 



(13) 
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Note again that we utilize that subset of the edge simplicity constraint corresponding to the averaged 4d-dihcdral 
angles. This simply halves the dimension of the phase space down to 3Nf. 

We also wish to have the symplcctic structure on the reduced phase space, which involves calculating the 
Dirac bracket: 

{•••}£ = {;■}- {;Si} A-'ij {£/,•}, (14) 

L L L L 

where Si denotes elements of the above constraint set and Ai.j— {Si, Sj} is the corresponding Dirac matrix. 
We present the details in Appendix C, and give a summary here. It emerges that for the basis parameters, the 
Dirac bracket is merely a rescaling of the corresponding Poisson structure: 



{A + ,A+} L 


= {A+,A + } 


= 


{A+,9+}l 


1 + l \A+ 
27 1 + 


h) = 


1 


{A+,4>+}l 


= {A+M 


= 


{4>+,0+}l 




8+} = 


f(i)({4>+}) 
A+ 


{</>+, <P+}l 


7 + 1 

2 {<P+,<f>+} 


sin a + 
= 7 A + 


{9+,9 + }l 




h) = 


l/(2)({0+}) 

7 A+ 



Remark: Let us reiterate that the commutation relations presented above are among the basis parameters 
and as such, they contain all the information about the symplcctic structure on the reduced phase space. One 
should not expect, however, that the Dirac bracket among more general quantities is a simple rescaling of its 
corresponding Poisson bracket. This stems from the fact that for the basis parameters the brackets are rescaled 
in different ways. To be a bit more explicit, let us calculate the Dirac bracket between a 3d-dihedral angle <j>+ 
and a 4d-dihedral angle 9+ (not the average, rather one with dependence on an edge of the shared triangle): 
{</>+, 6+}l- We know that 0+ = 9+ + ^ a+ and that every a+ may be written as a function of the 3d-dihedral 
angles. 5 Thus: 

= {0 + ,0+}l + {0+,/({0 + })}l 

= 2±i {0+ ^ +} + 2+l { ^ +)/({0+})} (i6) 
= 2±i { + ,0 + } + 2^_l w+ , /({0+})} . 

So the Dirac brackets involving non-basis 4d-dihedral angles are not simple rescalings of their Poisson bracket 
counterparts. Having said that, one finds that analogous commutation relations occur in the loop gravity phase 
space below. 



4.3.1 Relation to loop gravity phase space 

We claim now that we are at the phase space of loop gravity restricted to the graph topologically dual our 
triangulation A3. To cut to the chase, the subset of the loop gravity phase space is parameterized by a discrete 

5 To be a bit more precise, we are considering the situation: {<t>ijk+i®ij,kk' where we may write: 

Qij;kl+ = — a ij;kk' + a ij;kk" + a ji;Vl ~ a ji;l" I ■ 

s ' 

+/«>+» 
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triad and discrete su(2) connection: (xij r , rriij st ). With regards to the continuum theory, m.y is the holonomy 
formed from the Ashtckar-Barbero connection. The symplectic structure is then given by: 

{mij ru ,Xij s } L G = 7e" m,/ u , 
{Xif ,Xij S } LG = -fe rs t x ij t , (17) 
{m t3 ru ,m t] sv } LG = 0. 

As usual, we reduce to the su(2) gauge-invariant subspace. This is 3Nf dimensional and a basis is parameterized 
by: 





# in A 3 


# USED 


Areas: 


Aij 


N f 


N f 


3D-DIHEDRAL ANGLES: 


4>ijk± 


6Nt 


2N t = Nf 


Averaged 4d-dihedral angles: 


0ij± 


Nf 


N f 



These su(2) quantities are defined in an analogous fashion to their so (4) counterparts earlier, and we have 
utilized identical labels to facilitate comparison of the commutation relations (15) and (19). Upon calculating 
the commutation relations among basis parameters, one finds that: 



{A+,A+}lg 


= 


{A+, 


9+}lg 


= 7 




{A + ,<P + }lg 


= 


{>+, 


0+}lg 


/(i 

= 7 — 


({0+}) 
A+ 


{<!>+, <I>+}lg 


sina + 
= 7 A + 


{0+, 


9+}lg 


- 7 J ( 2 


A+ 



As one can see, these relations are quite similar to the ones we found in (15). There is a difference in their 
scaling w.r.t. 7, however, which can be explained in the following fashion. The LQG connection involves 7, 
schematically: A = T + jK, whereas we started with an SO(4) connection that was 7-indepcndcnt. As a 
result, the dihedral angle 8lg also involves 7. Indeed, replacing 9lg — > ^lg. °ne subsequently finds agreement 
between the brackets (15) and (19). 

Moreover, in the loop gravity phase space, it is rather the rescaled 4d-dihedral angles ^Qlg that are related 
via the 2d-dihedral angles ollg- Schematically: \Qlg = h®LG + Yl a LG- Thus, one finds that the commutation 
relations involving non-basis 4d-dihcdral angles are identical to those in (16). 

4.4 Stage 4: gluing constraints 

It is now time to deal with the remaining 3A^/ — 2A^ e edge simplicity constraints, which we refer to in their 
many forms as gluing constraints. It has been shown [52] that although the edge simplicity constraints relate 
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chiral sectors, they may been rewritten as relations among objects of one and the same sector. Imposing them 
is equivalent to imposing: 

%;fc/+ = %;fc'/'+ Or Otij-.kk'+ = Otji ; l>l+ ■ (20) 

In words, the left hand relation states that the different definitions for the 4d-dihedral angles agree, while the 
right hand one states that the definitions for the 2d-dihedral angles agree. The problem is that all the above 
variables are complicated functions of the basis parameters and thus rather tricky to impose, even via symplectic 
reduction. One may devise yet another equivalent set of constraints by the following argument. 

So far, the 3d-gcomctrical information has been encoded in both the area variables and the 3d-dihedral 
angles 4>ijk- However, for discrete 3d-gcomctries, the area variables provide already an over-complete set. Thus, 
the 3d-dihcdral angles on the geometric constraint hypersurface are actually determined as function of the areas: 

F ijk = (t> ijk + ~ $ijk+[{A}} . (21) 

The advantage of rewriting the constraints in this manner is that the Poisson (or Dirac) brackets between these 
constraints is straightforward to determine. Moreover, it leads to a Dirac matrix that factorizes over tetrahedra. 
In other words, these constraints are second class, while the two constraints per tetrahedron that fix the two 
independent dihedral angles form conjugate pairs. 

For the discrete spatial manifold described by the boundary of a 4-simplex, these are a necessary and suf- 
ficient set of constraints to reduce the phase space to the geometrical one. However, for regular triangulations 
that are larger than the boundary of the 4-simplex, there are generally more area variables than length vari- 
ables. Consequently, the constraints (21) are necessary but insufficient for such triangulations and one requires 
additional constraints that fix the areas as functions of the length variables. These latter constraints can be 
reconstructed using our basis variables if the full set of simplicity constraints holds [51]: 

A:A ij+ = aij [{l(A)}}. (22) 

Here 1(A) denotes some particular choice of how to construct the lengths from the areas A+. while aij(l) is the 
standard area of a triangle {ij} expressed as a function of its three edge-length variables. As we shall comment 
below, this set of constraints is 1st class, even Abelian, with respect to the Dirac brackets - both the brackets 
in (15) and the Dirac brackets after imposing the constraints (21). 

Thus, we now have an alternative set of 27V/ — 7V e constraints, of which a subset of cardinality Nf — N e is 
1st class, that capture the remaining 37V/ — 2N e edge simplicity constraints: 







# in A 3 


# USED 


Gluing: 


Fyfc 


= <t> ijk+ -$ ijk+ [{A}} 


Wt 


2N t = N f 


Area-length: 


A 


■ A ij+ = a ij+ [{l}} 


N f -N e 


Nf-N e 



4.4.1 Case 1: Boundary of a 4-simplex 

In this scenario, we are left with just the F constraints on the 3d-dihcdral angles. There are 27V/ = 10 such 
constraints, which rewrite the 10 basis 3d-dihcdral angles in terms of the areas. The bracket on the reduced 
phase space is: 

{;-}r = {;-}l-{;F} l A^ 1 {f, -} l , 
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where A 1 is the Dirac matrix corresponding to the constraint set F. With a little effort, sec Appendix D, one 
can show that: 



{A +1 A+} R = {A +1 A+} L = 



{A+A+}r = {A+,<t>+} L = 



{<t>+, <f>+} R = 



{A + ,9 + } R = {A + ,9+} L 



= 1 



{<P+,9 + }r = {*+[{A}],9+} R = 



OA 



{0+,O + }r = {0 +1 6 + }l-{0 + ,f}l A" 1 {f,0+}l 



We note that we have made significant progress on all the commutation relations except the final one containing 
two 4d-dihedral angles. One can make yet more progress via formal manipulations for this last bracket. However, 
ultimately one is forced to explicitly calculate a matrix inverse. Even for the simple case of a generic 4-simplex 
boundary, this is rather involved, see Appendix D. A full analytical evaluation is not possible, as this would 
require an analytical expression for the length of one simplex expressed as a function of its areas. Such a general 
expression is not available, as this would involve the general solution of a higher order polynomial. Having said 
that, we performed an explicit evaluation of the matrices in the equilateral case (see Appendices D for details) 
and found that it vanishes: 

{e + ,e+} R = o. (24) 

We can devise another argument to show that this bracket must vanish. On the initial SO (4) phase space, one 
may impose flatness constraints F e = id, demanding that holonomies around edges arc equal to the identity. 
Such constraints impose the dynamics of BF theory and are 1st class. Geometrically, the flatness constraints 
generate translations in the bi-vectors X, that is, they change the spatial (geometric) configurations. Then, let 
us reduce the initial SO(4) phase space to the one only allowing geometric configurations, that is, to the one in 
which we can consistently reconstruct the edge-length variables and dihedral angles. For the boundary of the 
4-simplex, we know that (4d) flatness can still be imposed and that the constraints take the form [60, 51]: 



cos%+ - cos&ij({A ij+ }) , 



(25) 



where 0ij({-A»j+}) is the geometric 4d-dihedral angle as computed from the areas of a flat 4-simplex. 6 The flow 
of the flatness constraints is now restricted to geometric configurations. Indeed, it leads to 10-dimcnsional gauge 
orbits. These correspond to the possible changes in the 10 edge-lengths, along with the associated changes in 
the 4d-dihedral angles, subject to keeping the 4d-geometry of the simplex flat. In order to produce a consistent 
flow, the constraints (25) need to be 1st class. This is the case if and only if the 4d-dihedral angles commute: 



{Cij,Cki}i 



sin I 



+ sin V M + 



d& l3 {A) dO kl (A) 



dA 



kl+ 



dAi 



{0ki+,6ij+}R 



(26) 



Here the equality sign ~ indicates an equation which holds on the (flatness) constraint hypersurface. The 
first two terms cancel each other, due to the Schlafli identity, which ensures that S a = ■ Aij+Qij(A) is a 
generating function for the dihedral angles, that is, SSa- = Ylij ^Aij + &ij(A). 



6 There are global and discrete ambiguities in reconstructing the geometry of a simplex from the areas. Here however, we need 
only local reconstructibility. 
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As a consequence, the existence of a consistent reduction of the flatness constraints to the geometric sector 
requires that the 4d-dihcdral angles commute. This argument can be extended to a class of triangulations, 
called stacked spheres, who for arbitrary edge length (satisfying appropriate triangle inequalities) are locally 
cmbcddable into flat space. Furthermore it also holds for the variables associated to the edges adjacent around a 
four-valent vertex, as the constraints for such four-valent vertices again agree with constraints imposing flatness 



4.4.2 Case 2: General triangulation 

For a more general triangulation, we have constraints among the areas, in addition to the constraints fixing the 
3d-dihedral angles as functions of the length variables. Their number can be easily determined as the difference 
between the numbers of triangles and edges in a given 3d triangulation. 

These constraints involve only the area variables, which commute with each other and are therefore 1st class. 
Thus, we can reduce with respect to these constraints by finding appropriate observables that are invariant under 
the flow of these constraints. The phase space will then be co-ordinatized by these observables and the induced 
symplectic structure leads to the same bracket as before {-,•}«• One class of observables is obviously given by 
the edge-length variables l e (A) associated to the edges e of the triangulation. 7 Nonetheless, one must describe 
how to compute these edge-length variables in terms of the areas - there will be different possibilities but by 
definition these will lead to the same result on the area constraint hypcrsurfacc. A second class of observables 
is given by: 



Here at (I) is the area of the triangle t expressed as a function of the edge-lengths, where these edge- lengths are 
determined by our area variables A + . 9 t + is the dihedral angle associated to the triangle t. As can be easily 
checked, the observables p e commute with the area constraints (22) and moreover, they are conjugate to the 
length variables l e with respect to the brackets {•, -}r. In the case that the dihedral angles to commute with 
each other, one finds that the momenta p e also commute: 



The first two terms cancel each other, as the area derivative factors are only non-zero if the triangles t, t' include 
both edges e, e', that is t = t' . 

Let us note that one finds the same kind of momenta p e in a canonical framework for Regge calculus [50] , 
in which the Regge action serves as the generating function for the momenta. There the momenta p e commute 
by construction. 

Thus, starting from an SO(4) phase space and implementing the Gaufi and simplicity constraints, we found 
so far a reduced phase space parametrized by length variables and conjugated momenta, which agree in their 
geometric interpretation with the canonical momenta found for Regge calculus. Also, the commutation relations 
between length variables and momenta agree with those derived from a Regge phase space. What is more difficult 
to show directly is that the momenta commute. However, we can also argue here in a somewhat similar fashion 
to the previous case: since we found that the geometric interpretation of the variables in our reduced phase 
space and matched those in the Regge phase space, then the dynamics should also be imposed in the same way. 
In the general case, the dynamics can be described by discrete time evolution steps that are in fact canonical 
transformations [50]. For the special situations where topological constraints insist on a flat dynamics, we 

7 Here, we shall change notation slightly to facilitate our argument: e and t refer to the edges and triangles of the triangulation, 
respectively. 



[50]. 




(27) 




(28) 
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argued that its consistent implementation required commuting momenta. Here, we conjecture that there is a 
unique symplcctic structure preserved by the canonical dynamics - in this case, it is the one of the Rcggc phase 
space, for which the momenta commute. 

5 Discussion 

Let us briefly recapitulate. We used a classical phase space approach to shed light on the role played by the 
Barbcro-Immirzi parameter in discrete quantum gravity. In the continuum description, this parameter appears 
in a canonical transformation, which turns the momenta conjugate to the triad variables into (commutative) 
connection variables. However, imposing primary and secondary simplicity constraints in the continuum, one 
can construct alternative (but non-commutative) connection variables [33]. It is argued in [33] that a path 
integral quantization based on these connections leads to a quantum theory, in which the Barbero-Immirzi 
parameter does not appear. 

We started from a discrete classical phase space. Introducing a discretization, we replace connection variables 
by holonomies. This is also done in the (continuum) loop quantum gravity quantization and is the manner in 
which the Barbcro-Immirzi parameter manifests itself in the quantum theory 

Indeed, after we impose a subset of the simplicity constraints, we find a phase space that can be mapped to 
the phase space of twisted geometries, which itself arises if we restrict Loop Quantum Gravity to a fixed graph. 
The Dirac brackets between basic geometric quantities are 7-depcndcnt and this leads to the appearance of 7 in 
the spectra of geometric operators. This latter fact shows that the quantum theories based on different values 
of the 7 parameter are not unitarily equivalent - even if we start from a finite dimensional phase space, as arises 
here from the restriction to a fixed graph. 

This first set of simplicity constraints imposes that for a given basis of gauge invariant phase space variables, 
the left handed-basis variables agree with the right-handed basis variables (for the 4d-dihedral angles up to 
sign). However, this does not mean that all geometrical variables constructible in the left-handed sector agree 
with their counterparts in the right-handed sector. It is rather after reducing with respect to a larger set of 
simplicity constraints that the phase space truly possesses the property: "left-handed geometry = right-handed 
geometry" . Moreover, after this further reduction the phase space dimension and Dirac brackets agree with 
their counterparts in the Regge geometries phase space and do not depend on 7. 

Thus, one has the choice to base quantization on either of these phase spaces. Due to the availability of 
a connection (or rather group holonomies), the loop quantum gravity-like phase space is given by a certain 
tensor product of T*(SU(2)) spaces and quantization is straightforward [61]. Moreover, the connection to the 
continuum theory via a projective limit is well developed [21]. In contrast, although a heuristic quantization of 
the Rcggc phase space is certainly possible, the imposition of triangle inequalities on this phase space is an open 
problem. On the other hand, a discrete (classical) canonical dynamics for the Regge phase space is available 
[50], whereas a dynamics for (discrete) classical twisted geometries has so far not been discussed. Indeed, the 
pertinent question here concerns what kind of dynamics one should impose on the non-geometric configurations, 
or in other words, whether these non-geometric configurations can be suppressed dynamically. 

Since the phase spaces have different sizes, one would expect different quantum theories, even at the dy- 
namical level. In particular, any appearance of 7 in the final dynamical predictions would be an effect due to 
the larger phase space of discrete loop quantum gravity, as compared to the Regge phase space. 
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A Conventions 

The discretization of 4-dimensional SO(4) BF theory relies heavily on the isomorphism between the space of 
bi- vectors (anti-symmetric second order tensors) and the algebra so (4). In our conventions, the commutation 
relations among the generators of the algebra take the form: 

[J a, Jb] =C ab c J c . (29) 

The indices are A = aa where a, a € {0, 1, 2, 3}. The structure constants are: 

C ABC ._ e Ar se Bs tS Ct r 

_ gab gabC _|_ gab gabC gab gabC gab gabC 

where the Kronecker delta and the summation convention on this space are given by: 

gAB ._ gabgab _ gab gab ^ 



(30) 



T AB U B C := X -T Ahl U hl c so that e AB e B C = S AC . 



(31) 



The Killing form is: 

C A CD C BCD = S AB (32) 
while the Hodge map e : so(4) — > so(4) so that: 

e(X) A = e A B X B (33) 
Using the Killing form and the Hodge map we can define two independent bilinear forms on so(4): 

< X, Y > := X A d AB Y b and -< X, Y y := X A e AB Y B . (34) 
We also define "scalar" and "bi-vector" products: 

X ■ Y := X A S AB Y B and {X x Y) A := C A BD e D c X B Y c . (35) 



We can use the Hodge operator to split the algebra into self-dual and anti-self-dual sub-algebrae: 

J± A ■■= P± AB Jb, where P ± AB := ^(S AB ± e AB ). (36) 

It is straightforward to check that these projectors are orthonormal, that is P AB P BC = 5 ss iP AC , furthermore 
P + AB + P_ AB = S AB . We may then proceed to generate all manner of projected quantities: 

X± A = P ± A A , X A ' M± AB = P ± A A , M A ' B = P ± B B , M AB ', 

C± ABC = P± A A , C A'BC = P± B B , C AB'C = p ± C Ql qABC' = p ± A ^ p ± B & p ± C Ql (jA! B' C 

(37) 

The following useful identity for the so(4) structure constants: 

C± ABC C± A'B' c = 2 3 { p ± AA' p± BB> _ p± AB> p± A' B) _ (3g) 

is in close analogy to the relation e abc e a b c — (8 aa S bb — S ab S a ) for the SO(3) structure constants. 

In our conventions, e A ' tl/p is a totally antisymmetric tensor density of weig ht 1 with e 0123 = 1. On the 
other hand, in absence of a metric, we define e\^ up as totally antisymmetric tensor density of weight —1 with 
£oi23 = !• This means that e XfJ -" p e\^ up = 4!. Furthermore, we define e IjA: := e 0ljA: and eyj. := eoijk- 
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B Gauge-invariant quantities 
B.l Definitions 

Let us first detail the gauge-invariant variables that occur in the main text and appendices: 



Areas: 



Squared volumes: 



Edge Lengths: 



2D- DIHEDRAL ANGLES: 



3D-DIHEDRAL ANGLES: 



4D-DIHEDRAL- ANGLES : 



.4 



ij± 



ijkl± 



Hjk± 



cos a %j -ki± 



COS (pijk± 



cos tj - k i± 



1 

71 



Xi i± \ 



[Xij ■ (Xik x Xii)\ ± 
\N ijk \ 



i 



j ± 



Wij k \\N l3l \ 
■ Xik 



2 Aij Aik 



Nijk ■ MijNju 
\Nijk\\N ja \ 



sm aij. k i± := 



sm <p ijk ± := 



sm Oij-ki± 



±Ai,Vi 



ij 'ijkl 



\NvkWNw 
2 2 Aik 



J i 



AAjj(X ik - M t] Nju) 
\Nijk\\N ja \ 



where N ijk ± := (Xy x X ik )±. 

Remark: We define both the cosine and sine of each angle above. Since they both involve square roots, we do so 
in order to specify their signs. A priori, this is an ambiguity, although it is removed by picking one particular sign 
consistently. Moreover, for sin^±, we wish to reflect the fact that 9 + = —9- for geometric configurations. Thus, 
we pick the factor X ik ± ■ Mij±Njn± in the numerator, which satisfies X ik+ ■ Mij + Njn + = —X ik ^ ■ Mij—Nju- 
on the constraint hypersurface. Also, sin a is essentially a tri-linear function on phase space and the presence 
of V (rather than |V|) represents the change of sign of smctij-ki as one permutes the indices k and I. 

Remark: Let us also details the evaluation of these quantities for geometric configurations. In this situation, 
there is a reference frames attached to each tetrahedron. Pick a vertex of the tetrahedron as the origin and use 
the three edge-vectors ei.2,3 emanating from it as the basis. Taking the discrete bi- vectors to be: 



Xij A 



Xa A 



.Abb , b „ b 

e ei e2 , 



e Abb e 3 b ei b 



Lift 



.Abb „ b „ b 

e e2 e3 , 



(39) 



v A —Y A _ y A _ y A 

^ira ^ij -™-ik j 



leads to the following: 



Aij + 
Vijki+ 



1 



ei\\e 2 \ sm a ijt ki+ , 



2 

18kl 2 



v, a 



^~ .abed „ b - c ~ d 

—e ei ei 
6 



(40) 



Note that w ? a satisfies: \vA = volume of tetrahedron. 
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B.2 Providing a basis 

Given a triangulation A3, we claim that we can parameterize the gauge-invariant phase space using, in each 
chiral sector, all the areas, two 3d-dihedral angles per tetrahedron and one 4d-dihedral angle per triangle. These 
is a set with 6Nf elements, which is certainly the correct number but let us now show that we can reconstruct 
the rest of the gauge-invariant quantities from these basic set. 

Let us first deal with intrinsic quantities, that is, quantities pertaining to a single tetrahedron marked {i}, 
with four faces {ij}, {ik}, {il} and {im}. From our proposed basis, we have the following parameters: 

Aij , Aik , An ; Aj m 1 cos (frijk , cos (\>%ji . (41) 

We suppress the ± here for ease of reading. Then, the other 3d-dihedral angles are given by: 



1 Aij 
COS (pijm = -—. < 

1 / A-'ij A-ifc A-ai Aim, A-ijAn 

cos0 ifem = - — -j — + — — + 2cos^i - 



Aik 




An 


hjk A 


- COS fajl 


A ■ ' 


Aik 


J2 


A ■ 

n ™ 1 


A ■ 


AikAi m 


A ik 


A 2 

A ik 


An 


A - 


A; A- 


A 


An 



2 \ AikAi m Aim AikAi m Aik AikAim 



(42) 



, 1 1 Aij Af k An Ai m , n 1 AijAik 
C0S<t>llm = 2 \ T7a~ + ^~^^~ — + 2cos^ fc 

, 1 J ^^ij Aik An Aim ^ ^ ^jy ^ Aij ^ i _ Aij j 

2 I Aik An An Aik Aik An 3 An 3 Aik J 

Moreover, the 2d-dihedral angles may be constructed from these 3d-dihcdral angles: 

COS (f>ikl - COS (f>ijk COS 4>iji . , 

cos aij-ki = ^— 7^— — • (43) 

sin (pijk sin (piji 

We are left now with extrinsic quantities, that is, the 4d-dihedral angles. Consider a triangle {ij}. It is 
a face of two tetrahedra {i} and {j}. They have faces {ij},{ik}, {il}, {im} and {ji}, {jk 1 }, {jl'}, {jm 1 }, 
respectively. Thus, for the triangle {ij}, there are three 4d-dihedral angles: 9ij ; kk', Oij;W and 6ij- m m'- Using 
the relations: 

Oij-w = Oij-kk 1 + oiij-ki + otji-k'i' , 

(44) 
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B.3 Commutation relations 



We provided a very cursory digest of the commutation relations among the gauge- invariant parameters of interest 
in the main text. We provide a more complete list here: 



{Aij± , A mn ±\ — , 

{Aij±, 4>mnp±\ = , 

7 1 

{4>ijk±,4>ijm±} = ±2 — — Smaij-km± , 

7 ± 1 Ai j± 



7± 1 ' 

7 / esc <fijjk± cot 0jjfc± 
7 ± 1 I A lk ± A lJ± 



7 1 

{4>ijk±, 9ij;mn±} = ±2 — — COt <j>ij m ± COS Otij;km± , 

7 ± 1 Aij± 



7 1 

{fe±,L;jiii) = T2 — CSC (j)imj± COS Oiij- } km± 

7 ± 1 Aj 3 -± 



7 1 / 

{^;iIi,C.j;mni} = ±2 — — COt <f>ijk± COt 0y m ± SU1 O^jfcmi + COt <^± COt <f>ji n ± sin aji ; / n ± 

7 ± 1 Ay± V 



7 1 

{%;fei±,^m;i'n±} = T2 — -j COt (f)ijk± CSC <f>i m j± sin (Xij- } mk± 

7 ± 1 Ay-t 



7 1 

{dij;kl±, 6im;kn±} = T2 — CSC 0ijfe± CSC 4>imk± Sm a l k;jm± 

7 ± 1 Aife± 



{9ij;kl±: Oik;jn±} — 
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C The calculation of {-, - } L in STAGE 3 



As we mentioned in the main text, we intend to explicitly calculate, via symplectic reduction, the Dirac bracket 

L 

on the phase space reduced by the constraint set Si= {d, C, e}. At first glance, this might appear to be a rather 
arduous task. The exists a neat fact, however, that simplifies the calculations immensely. Since the constraints 
are simple linear equalities between pairs of phase space parameters, there is a remarkable similarity between 
the commutator of constraints and commutators of phase space parameters, as illustrated in (45): 



A = {d,e} + 
B = {C,C} + 
C = {c,e}+ 
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7-1 

2 

2-) 



{A+,9+}, 

{<£+,</>+} , 



(45) 



7' 



v 



{e,e}_ 
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where {•,•}+ denotes the Poisson bracket evaluated on the reduced phase space. As a result, one may derive 
the symplectic structure on the reduced phase space from the block form of the inverse of the Dirac matrix. 
The Dirac bracket is defined in (14) as: 



{vH = {•>•}-{•, Si} K^ij {Sj,-} , 



where: 



L 



/ A \ 



B C 
\ -A -C T V J 
As a simple example, let us calculate: 



L 1 



( A- 1 £ T B- 1 £A- l +A- 1 gA- 1 -A- 1 £ T B- 1 -A~ l \ 



-B- X £A- X 
A- 1 



B- 1 









(46) 



{A.. + ,9..+} L = {A. + ,0..+}-{A. + ,e}(A- 1 ) 3 i{d,0.. + } 



T-U-C^l A A' 1 A 



2l 

7 2 -l 
4 7 2 

7 + 1 
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^4 



27 



{A.. + ,6..+} 



(47) 



The other brackets follow in a similar fashion with the results summarized in (15). 
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D The calculation of {-,-}r in STAGE 4 

We shall reduce with respect to the constraints {f, a} for the case of the boundary of a 4-simplex. 



D.l 4-simplex boundary 

For the case of a 4-simplex boundary, the remaining constraints simplify to the 2nd class set {f}, and so the 

R 

Dirac matrix simplifies to A = {0+, 0+}- Thus, its inverse is easily calculable as the matrix elements satisfy: 



R 



1 



R 

Aji 



(48) 



When explicitly calculating large Dirac matrices, it becomes rather important to choose the constraint set wisely, 
as this will cut down the amount of work considerably. An independent subset of constraints was found in [?] , 
which we shall utilize here. First of all, let us choose a parameterization of our 30-dimcnsional gauge-invariant 
phase space: 8 



A12+, Ai 3+ , A23+, A 2 4+, A 3 4 + , A 35+ , A i5+ , A 4 i + , A 5 i + , A 52 +, 

0125+ , 0135+ , ^231+ , 0241+ , 0342+, 0352+, 0453+, 0413+ , 0514+ , 0524+, 
#12;55+, #13;55+, #23;11+, #24;11+, #34:22+, #35:22+, #45;33+, #41;33+, #51;44+, #52;44+, 

Then, an appropriate constraint set is: 

F 125, F 135, F 231, F241, F342, F352, F453, F413, F514, F524, 

Due to the nature of the relation (48), the inverse of the Dirac matrix is easy to compute: 
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where: 



A x 
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7Sinai5,23+ 



8 We have chosen here the 9ij,kl rather than the averaged angles 9jj for purely historical reasons. 
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(51) 
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Moreover, in one shifts all the subscripts by i (mod 5) in the previous equation. Then, we have: 
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where: 
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The final bracket to calculate is: 



{6+,0+U 






(1) S1235 


(2) 

92315 


f2) 

92415 














(2) 

~ V '31254 


- (2) 92154 


- (1) S1235 





(2) 

93215 








(2) 

— 93152 





(2) 

91435 


- <2) 91354 





- (2) S2315 


- (2 'S3215 





(1) 92341 


<2) 93421 


<2> 93521 














- (2) S2415 





- (1) 92341 


(1 


<2) 94321 








- (2, 94213 





(2) 

V '92541 








- (2) 93421 


(21 

- 94321 





93452 


(2) 

94532 


(2) 94132 











(2) 93152 


- (2) 93521 





- (1) 93452 





(2> 95432 








— (2) 95324 














- (2) 94532 


- <2) 95432 





(1) 94513 


<2) 95143 


<2> 95243 





(2) 

— 31435 





f2) 

V ' 94213 


" (2) 94132 





- (1) 94513 





<2 >91543 





(2 ' 91254 


(2) 91354 














(2) 

— 95143 


- (2) 91543 





(1) 95124 


(2 'S2154 








- (2) 92541 





<2> 95324 


- <2> 95243 





- (1) 95124 






where: 



(i) 



1 1 



gijki 
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1 A lk+ 
1 1 



9ijkl 



csc (pijk+ esc 4>uk+ sin a ik -ji+ , 
cot (j)ijk+ csc 4*uj + sin a ij: ik+ ■ 



7 A ij+ 

In choosing the constraints as we did, we have postponed the pain to the following evaluation: 

d$+(A) 

{f,Ml = {<t>+,o + } L - -^-^{A + ,e + } L , 
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While we cannot specify explicitly, we can in principle circumnavigate this obstacle for the 4-simplex 

boundary since the 10 areas actually already specify a set of 10 edge lengths {/ + } for the boundary in 

question. Let us remark that these lengths {1+} are not identical to the {L + } defined in Appendix B.l at this 
stage, although they do coincide on the geometrical phase space. Moreover, the may be written as explicit 
functions of these edge lengths, that is, as $ + (A(7)). This manipulation allows us to calculate: 



d$+(A) d$+(A(l)) dl+ d$+(A{l)) f dA_\ 



dA+ 



dh 



dA A 



dh 



(57) 



To clarify, for a 4-simplex boundary: 
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(58) 



hjkhjl 



where the hats record the fact that in the main text we were dealing with exterior dihedral angles, while here 
we are defining interior dihedral angles. Thus, $ = tt — <f> and a = n — a. Under this transformation, there 
are a number of important sign changes. Additionally, the subscripts ijklm = er(12345) where a e S5. The 
relevant derivatives are straightforward to calculate: 
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1 COt Otik-jm ~ TTJ— cos ®ijk COt a ij;mk , 



(59) 
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1 7 2 7 

1 lijkHlm 



4 Mi Mk 

While it may seem plain sailing from now on, one should remember that we have yet to invert the matrix 
While some progress can be made for general configurations, we must admit that we have not found a way to 
simplify the sums of products of cotangents. Thus, we specialize yet more to the equilateral configuration in 
order to perform an explicit inversion. 

D.2 The equilateral 4-simplex 

8A(l) dl(A) 
Now, we must invert the matrix — — — to get the matrix of partial derivatives: n - . We shall order the areas 

dl OA 

and lengths as: 

A 12 A 13 A 23 A 24 A 34 A 35 A 45 A 41 A 51 A 52 , 

(60) 

ll25 ll35 ^231 ^241 ^342 ^352 ^453 ^413 ^514 ^524- 
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Furthermore, we shall consider only the equilateral case in which: 
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Then, the matrix evaluations yield: 
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while: 9 
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Ultimately, one finds: 

9+}r 

which is the result we were expecting. 
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